The many-body Green function of degenerate systems 
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A rigorous non perturbative adiabatic approximation of the evolution operator in the many-body 
physics of degenerate systems is derived. This approximation is used to solve the long-standing 
problem of the choice of the initial states of Ho leading to eigenstates of Hq 4- V for degenerate 
systems. These initial states are eigenstates of PqV^Po, where Po is the projection onto a degenerate 
eigenspace of Hq. This result is used to give the proper definition of the Green function, the statistical 
Green function and the non-equilibrium Green function of degenerate systems. The convergence of 
these Green functions is established. 
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Non-perturbative Green function methods, such as 
the GW approximation [1] or the Bethe-Salpeter equa- 
tion [3,13], have brought remarkable progress in the calcu- 
lation of the electronic structure and dielectric response 
of semiconductors. The extension of these methods to 
transition metals systems faces a serious difficulty: the 
standard Green function can only be defined when the 
initial state |0) of the system without interaction is a 
single Slater determinant. In physical terms, each single- 
particle orbital or Bloch state has to be either occupied or 
unoccupied at zero temperature. However, the physics of 
transition metals often contradicts this requirement. For 
the example of a V'^+ ion in an octahedral environment, 
we do not know a priori how the two Sd electrons are 
distributed over the six degenerate t2g orbitals (with up 
and down spins). 

More generally, for a system described by a Hamil- 
tonian H = Hq + V where the ground state of Hq is 
degenerate, we need to determine the parent states, i.e. 
the initial states of Hq that evolve into eigenstates of H 
by adiabatically switching the interaction. 

Degenerate systems being ubiquitous in quantum 
physics, this long-standing problem has been discussed 
in chemical physics [1, 0] , nuclear ph ysic s @, 0| , atomic 
physics JSk |9|] and solid state physics [lO| . Esterling and 
Lange [lO| summarized the situation as follows: "Since 
Hq has degenerate ground states, the choice of the state 
|0) must be made with care, and this may be considered 
the key to the problem." This question is also crucial 
in many-body physics because the Green function of a 
degenerate system has to be defined from a parent state. 

In the present paper, we give a simple method to ex- 
plicitly determine the parent states and to define the 
Green function of degenerate systems. Through a non- 
perturbative analysis of the evolution operator of a de- 
generate system, we determine the exact form of its sin- 
gularities. This enables us to derive: (i) an easy and 



explicit method to determine the parent states; (ii) a non- 
perturbative proof that the Gell-Mann and Low formula 
generally converges only for these parent states; (iii) the 
formula for the Green function of degenerate systems; 
(iv) the validity of the so-called statistical Green func- 
tion; (v) the singularity structure of the non-equilibrium 
Green function of degenerate systems. 
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Adiabatic switching. Many-body theory 
usually based on the adiabatic switching of the inter- 
action, i.e. the transformation of the time-independent 
Hamiltonian H = Hq + V into the time-dependent one 
Hq + e~'^^*'^V. Adiabatic switching turns the non degen- 
erate ground state |0) of Hq into an eigenstate |\E'gml) 
of H first proposed by Gell-Mann and Low flS^] in 1951 
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where the evolution operator U^{t,t') is the solution of 



.dUe{t,t') 



s*^°*e--'l*lye-*^»*[/e(t,i'). 



with the initial condition U^{t' ,t') — 1. The wavefunc- 
tion Iv^gml) is then used to build the Green function of 
the system 11, [l3|. However, Gell-Mann and Low did 



not prove that the limit of eq. ([T]) exists I4I. The con- 
vergence of I^gml) for nondegenerate systems was first 
established by Nenciu and Rasche in 1989 Ii*!. 

For a degenerate ground state |0) of Hq, the Gell-Mann 
and Low formula generally fails to converge when e — > 0, 
as can be seen even for a trivial two- level system [l^ . In 
the following, we use recent advances in the mathematical 
analysis of the adiabatic approximation (see [l6j for a 
review) to extend the Gell-Mann and Low formula to 
degenerate systems. 

Adiabatic approximation. In this section, we set up the 
notation and give the theorem that enables us to calcu- 
late the parent states and to define the Green functions. 
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FIG. 1: Exemple of allowed eigenvalue pattern 



We consider i < and we rewrite the time-dependent 
Hamiltonian as H{s) — Hq + e'^V, where s = et is the so- 
called slow variable. The eigenvalues of H(s) are denoted 
by Ej{s) and its eigenprojectors by Pj{s). We recall 
that, if the eigenvalue Ej{s) is rij-iold degenerate, the 
eigenprojector is Pj{s) = YlkLi\'P]k{s)){(pjk{s)\, where 
il'Pjkis))} is a set of nj orthonormal eigenstates of H{s) 
for the eigenvalue Ej{s). For notational convenience, we 
denote Pj{—oo) by P^ in the rest of the paper. The model 
space M is the vector space generated by the eigenstates 
corresponding to A'o eigenvalues of Hq = H{—oo). Each 
eigenvalue of Hq can be degenerate and is possibly split 
by the perturbation V, so that the Nq eigenvalues of Hq 
become N eigenvalues Ei{s), . . . , E]y{s) of H{s), with 
N > Nq. Each Ej{s) can be degenerate and the eigen- 
values are allowed to cross (fig. [Ij. For an octahedral 
ion, we have A'o = 1 with degeneracy 15, and there 
are = 4 interacting states: ^Aig, ^Eg, ^T2g and '^Tig, 
with degeneracy nj = l, 2, 3 and 9, respectively. 

A key tool of our approach is A(s,so), the rotating 
frame operator [13, [3 , that relates the eigenstates at sq 
and s: A{s, so)\(f>jkiso)) = |0ifc(s)), so that 



A{s,sq)Pj{so) = Pj{s)A{s,sq) 



(2) 



Using standard technical assumptions [l9|, we recently 
obtained [1^ a rigorous approximation of the evolution 
operator projected on each eigenspace: 



f/e(0, -oo)P" ~ e^^^/-'A(0, -oo)FO 



(3) 



where 9j = — J_^iEj(s) — Ej(—oo))ds. In particular, 
the divergences of the evolution operator are entirely de- 
scribed by the factor e*^^/^. 

Construction of the parent states. The parents states 
are the eigenstates \(f)) of Hq such that [4(0, — oo)|0) 
tends to an eigenstate Ivt) of H, up to a (divergent) 
phase. Therefore, the parent states are naturally de- 
fined in terms of C4(— cxd, 0)|5') and it seems that the 
interacting states l^*) are needed to define the parent 
states [1, 1^ im . We now show that the parent states 
have a more simple and explicit definition as eigenstates 
of Pj and we explain how can be calculated by stan- 
dard time-independent perturbation theory. 



For notational convenience, we denote e** by A and the 
eigenvalues and eigenprojections are written in terms of 
A. We denote by Ej{X) and |^jfe(A)) the eigenvalues and 
eigenstates of Hq + XV (so that Ej{X) — Ej{s)). They 
can be expanded as |20j 
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with the normalization {'p'jklVjkW) = 1- The eigen- 
states |^j7j(A)} are assumed orthonormal only at A = 0, 

where P," = EEil^°.)(^°J- ^ , 

The time-independent Schrodinger equation 

{Ho + XV)\^,k{X)) = E,iX)\^,,{X)), 

gives, to order 0, Hq\lp^-,,) = E'^\ip^^), so that \ip^f.) is 
an eigenstate of Hq with energy E^. We assume that 
is one of the Nq eigenvalues of the model space, 
so that I'/Sjj,) belongs to the model space. However, 
the degeneracy of Ej as an eigenvalue of Hq is gener- 
ally larger than the degeneracy of Ej{X) and we need 
more information to determine the Uj states |<Pjfc). The 
Schrodinger equation to order A gives us {HQ—Ej)\(pj,^) = 
{Ej — V)\(p^i^). This equation can only have a solution if 
«|(4 - V)\ip%) = 0, where is a complete set 

of eigenstates of Hq with energy Therefore, the ini- 
tial states Iv^^fc) are eigenstates of Hq with energy Ej and 
eigenstates of Pe-oVPeo with eigenvalue E}, where Pro 
is the projection onto the eigenspace of Hq with eigen- 
value Ej. In general, the degeneracy is split at this or- 
der, in the sense that there are only nj states that are 
simultaneously eigenstates of Hq with energy E^ and of 
Pc-oVPc-o with eigenvalue E]. Otherwise, for instance 
when Ppjo T^Pro is zero by symmetry, the equations com- 

3 j 

ing from higher powers of A must be taken into account 
to determine |t/?°i.). In that case, the second order is usu- 
ally enough gJl^but methods have been developed to 
treat any order |23[ . 

We generally have no a priori knowledge of nj and Ej . 
However, we can calculate all the eigenstates of Hq and, 
for each energy E^, we can diagonaHze P^oVP^o. Then, 

3 3 

each state must be examined to see if it cannot be fur- 
ther split by higher order terms. When degeneracy is due 
to the symmetry of the Hamiltonian H(s), this can be 
deduced from the dimension of the irreducible represen- 
tations to which the states belong. The computational 
effort required to construct |<p^j.) is small because it is an 
eigenvalue problem in a vector space whose dimension 
is the degeneracy of Ej, which is small in applications. 
From the states \ we build the projector P? and we 
define a parent state as a state 10) such that, for some j, 



(4) 



In practice, the parent state is one of the \tf'^- 
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Generalized Gell-Mann and Low wavef unction. We 
show that the parent states previously defined lead to 
convergent Gell-Mann and Low wavefunctions. For a par- 
ent state \4>j) such that Pf\(f)j) = \<j>j), eq. ^ enables us 
to write 



f/,(0,-cx3)|^,) 



= C/e(0,-oo)P°|(/.,) 
~ e*''^/M(0,-cx))P°| 
~ e**^/M(0,-oo)|(/)j) 



Therefore, the following limit exists: 

f/,(0,-oo)|(/.,) 
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o(^,|C/,(O,-oo)|0,) 
A(O,-(X3)|0,) 
(0,1^(0,-00) 10,) • 



(5) 



The Gell-Mann and Low wavefunction I'I'gml) is in- 
deed an eigenstate H^ + V with energy i?j(0) because 
Pj(0)|*gml) = I^'gml)- To show this, we use eq. ([2]): 

P,(O)A(O,-(X3)|0,) = ^(O,-w)Pj'|0,) 
= .4(O,-«))|0j). 

In practice, we are interested in the Gell-Mann and 
Low wavefunction that is the ground state of Hq + V . 
How should we choose the initial state \(f)j) for this to 
happen? In the non degenerate case, it is often assumed 
that the ground state of Hq leads to the ground state 
of Hq -I- V . When degeneracy is due to the presence of 
symmetry, band crossing can occur and one should try 
the corresponding to the lowest energy E'^ for each 
irreducible representation. A typical example of band 
crossing in the presence of symmetry is given by Tanabe- 
Sugano diagrams of the multiplet theory [24|. For a 
small crystal field, the ground state has the highest spin 
value (Hund's rule), but as the crystal field parameter 
increases, a low spin state can become the ground state. 

Green functions. The expression for the Green func- 
tion is usually derived under the assumption that the 
ground state of is non degenerate [H, 12|. Our results 
enable us to determine how this expression is extended to 
the case of degenerate systems. Now we formally extend 
our previous results to Fock space. 

To follow the usual argument [13], we repeat the cal- 
culation by starting from a positively infinite time. In 
terms of the slow variable s = —£\t\, the switching func- 
tion e~^l*l is the same for positive and negative times. 
As a result, the rotating frame operator is the same 
but the divergent phase changes sign: [/^(O, -|-oo)|0j) ~ 
g-ie,/e^(0, -oo)|0j). Therefore, 



lim 



i7e(O,+oo)|0j) 



A(O,-(X3)|0,) 



(0j|C/,(O,+^)|0j) (0,|A(O,-(»)|0,)- 

In other words, the Gell-Mann and Low wavefunctions 
obtained from positive and negative infinite times are 



equal. This non-trivial result is due to the fact that the 
switching function f{t) = e"*^'*' is even. 

The two-point Green function is defined by 



Gix,y) = 



(*GML|Off|*GML) 



gmlI^gml 



where x = (r,i), y = (r',t'). Oh = T{i^H{x)i^Uy)) 
is the time-ordered product of fields in the Heisenberg 
picture and I^'gml) is defined by eq. ([5|). Standard ma- 
nipulations [12,] transform it into 



G{x,y) 



lim 



{(j)j\Ue{+CO,-Oo)\(l)j)' 



(6) 



where = U^i+oo, t)^{x)Us{t, t')ip^ {y)Ue{t' , ~oo) and 
= -U^{+oo,t')'4!'^ {y)Ue{t' ,t)ip{x)Ue{t, -oo) ii t > t' 
and t < t', respectively. 

The expression for the Green function generally con- 
verges only when the initial state is a parent state. In- 
deed, consider a state |0) in the model space and write 
it as 10) = \<f>j). where |0,) - pO|0). Thus, 



C/,(O,±«))|0)~^. 



A(O,-oo)|0,) 



If there is more than one j in the sum, the phases 9j are 
generally different (in the absence of eigenvalue crossing, 
they can be shown to be different). Therefore, the phase 
factors in the numerator and denominator of eq. ([6]) do 
not cancel and the expression has no limit for e ^ 0. 

Statistical Green function. The Green function of the 
previous section has a non-ambiguous meaning when \4>j) 
is the parent state of a non-degenerate interacting state. 
However, when the interacting state itself is degenerate, 
there is no reason to choose any particular parent state. 
To solve that problem, Layzer [25| defined the statisti- 
cal Green function as an equal-weight average over the 
degenerate states. Such a statistical Green function was 
advocated, for instance, by Alon and Cederbaum [j^ . 
The statistical Green function can preserve the symme- 
try of the system: in the example of a spherically sym- 
metric Hamiltonian, the Green function obtained from 
any state \£m) with i gives non spherically symmet- 
ric charge density, whereas the statistical Green function 
obtained from the mixed state \£m)(2£ + l)~^{£m\ 
gives a spherical charge density [23|. We are now able 
to prove that the statistical Green function has a well- 
defined limit when e ^ 0. 

To define the statistical Green function of a degenerate 
interacting system with energy -E'j(O), we use the den- 
sity matrix p — (l/nj) J^^Li \ f%){'P%\7 where the states 
|V5°fc) are those used to calculate Pj'. We assume that the 
degeneracy of Ej{0) is Uj. Then, 



G{x,y) 



lim 



Tr(pXe 



*o Tr p{/,(- 



-«))) 



(7) 
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If we put = in eq. ([6|), the divergences of the 

numerator and denominator are e^'^^/^, which does not 
depend on k. Therefore, the divergent phases of the nu- 
merator and denominator of eq. ([7]) are equal, and the 
statistical Green function is well defined. For our octa- 
hedral V^^ ion, the Green function is defined with the 
density matrix built from the nine degenerate states with 
^Tig symmetry if the ion is high spin. 

Non-equilibrium Green function. In the study of non- 
equilibrium systems, it is often convenient to run the evo- 
lution operator over a closed time path instead of taking 
the limit t +00. In this so-called Keldysh approach, 
the Green function G(x, y) is calculated by a formula in- 
volving no denominator [28|: it is the limit for e ^ 
of = (0|f/e(-oo,O)O//f7e(O,-oo)|0). As for the stan- 
dard Green function, this expression generally converges 
for degenerate systems only when is a parent state. 
To see this, we expand again a state of the model space 
over parent states: |0) = ^ ■ Then, 



'''•)/^(,/),|A(-^,0)O^^^(0,-oo)|</.,). 



This expression converges for e — > when there is a single 
phase, i.e. when is a parent state. Otherwise, the 
limit generally does not exist. 

Conclusion. The determination of the parent states 
and the proof of convergence break the last deadlocks in 
the determination of the Green function of degenerate 
systems. The main difference with the non-degenerate 
case is the fact that, for many degenerate systems, the 
parent state is not a single Slater determinant. To see 
this, consider a Hamiltonian where Hq is the restricted 
Hartree-Fock Hamiltonian of a 3d" transition metal ion 
and V is the sum of the remaining atomic Coulomb inter- 
action and of an effective potential representing the influ- 
ence of the surroimding atoms. Then, the parent states 
are exactly the cigenstates of the crystal-field Hamilto- 
nian and they are generally not single Slater determi- 
nants. In that case, the structure of the Green function 
is more complex because of the so-called initial correla- 
tions [29] coming from the matrix elements between the 
different Slater determinants. Perturbative [2^ and non- 
perturbative [lOl methods have been developed to tackle 
initial correlations. Finally, our approach also gives a 
non-perturbative proof of the convergence of the effec- 
tive Hamiltonian ^3L] . 
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